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Abstract 

In this paper we report the results of electrical resistivity (1.5< T < 300K) 
and point contact spectroscopy (PCS) measurements on single crystals of 
metallic sodium tungsten bronze with varying sodium content. We have shown 
that the electron- phonon coupling function as measured through PCS can ex- 
plain quantitatively the large temperature dependence of resistivity p seen in 
these materials over the entire temperature range. The electron-phonon cou- 
pling function shows predominately large peaks for phonon frequency range 
of 30 meV < cu < 100 meV which match well with the calculated optical 
phonons for WOg octahedron. The integrated electron-phonon coupling con- 
stant A from this data is ~ 0.25-0.45, depending on the Na content. 
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I. INTRODUCTION 



Electrical conduction in oxides, in particular, in metallic transition metal oxides is a topic 
of current interest0i. Over the last decade three classes of phenomena (namely superconduc- 
tivity, colossal magnetoresistance and insulator-metal transition) have kept the oxides in the 
center stage. For the high T c cuprates even the normal state resistivity with a linear tem- 
perature dependence is an unresolved issue. While the problem is most severe in the high T c 
cuprates the temperature dependence of the resistivity in the other transition metal oxides 
have a number of interesting features in store. To be specific in this paper we are concerned 
with normal (i.e, non-superconducting) metallic oxides which show the p(T) as in a metal 
(i.e, dp/dT > 0) and would like to investigate whether the temperature dependence of p can 
be understood quantitatively. In this context two important observations are noteworthy. 
First, in these oxides the resistivity even in the metallic state is quite high and second, 
the magnitude of the temperature dependence of resistivity, Ap(T)(= p(T)-p , where p is 
the residual resistivity) is also very high. In fact, incomparison to conventional metals and 
alloys, where the Ap (often <100/if2cm) arises from electron-acoustic phonon scattering, the 
Ap in the metallic oxides can be considered giant. If the source of the extra Ap in oxides 
is scattering of electrons by phonons then it is clear that there must exist a source of large 
density of phonons in these oxides. Such a contribution has been found in certain oxides. 
It was shown in oxides like Bai-^K^BiOa (superconductor with T c ^17K at x ~0.47) and 
Ndi.85Ceo.i5Cu04 (superconductor with T c « 22K) that the scattering from optical phonons 
in these materials make a significant contribution to the temperature dependence of p in the 
normal state!. These studies obtained the phonon data from the superconducting tunnel- 
ing spectroscopy which gives the Eliashberg function a 2 F{uj) and used it to calculate the 
resistivity. This investigation is a very good beginning in the quantitative understanding of 
the temperature dependence of the resistivity of metallic oxides. This study, however, had 
one big drawback. Most of the data were on polycrystalline thin film or bulk pellets. Also 
the samples used for the tunneling data were not the same as those used for the resistivity 
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studies. In order to improve upon this, we did a complete investigation with single crys- 
talline Sodium Tungsten Bronze system (Na^WOs) and used the same sample for both the 
resistivity measurement and the phonon spectra measurement from which the Eliashberg 
function a 2 F(u) can be obtained for calculation of the resistivity In this investigation we 
take a fresh look into the problem in order to see that the large temperature dependence 
of the resistivity seen in the metallic oxides indeed arise from the electron-optical phonon 
scattering. By making the investigation on a different and somewhat simple metallic oxide 
system we also establish the generality of the phenomena. 

Our experiment contains essentially two distinct parts: (1) measurement of the resistivity 
p over the temperature range 1.5 < T < 300K and (2) measurement of the phonon spectra 
(i.e, a 2 F{ui)) using the technique of point contact spectroscopy done at 4.2K. We then use 
the experimental phonon spectra to calculate the temperature dependence of the resistivity. 

The point contact spectroscopy (PCS) was first used by Yansoni to obtain the phonon 
spectrum (PS) of conventional metals and dilute alloys from the non-linear current volt- 
age characteristics of the point contact junction. This spectroscopy has the main merit 
that, it can give the function a 2 F{ui) which is simply related to the Eliashberg function 
a 2 F(u) (electron-phonon interaction function) for any metal including those with very weak 
electron-phonon coupling. The relation between a 2 c F(uj) and a 2 F(u) is essentially a factor 
which depends on the geometry of the junction formed between the tip and the sample!'! as 
discussed later on. This method enables one to derive the phonon spectrum of those metals 
which are not superconducting. We have thus used this technique to obtain the a 2 c F{oj) 
and use this function to fit the resistivity of the metallic sodium tungsten bronzes. It must 
be pointed out that this is the first time the phonon spectra is being studied using PCS 
technique in this type of metallic oxides. 



3 



II. THE SODIUM TUNGSTEN BRONZE SYSTEM 



The sodium tungsten bronze system, (Na x W0 3 ), investigated in this work is a 3- 
dimensional cubic perovskite oxide. The corner sharing WOg octahedra forms a 3- 
dimensional network in WO3 which is an insulator because the conduction band formed 
by the overlap of oxygen 2p and W 5d orbitals is empty. On addition of Na, it contributes 
one electron per Na atom and the conduction band gets filled. For x > 0.2 the system 
makes a transition to a metallic state. However, the cubic structure is stabilized only for 
x > 0.4. This material is very similar to Re03 which may be the most metallic of all oxides 
and in single crystalline form can be comparable to copper!. In our system for x = 0.9 the 
residual resistivity, po — 3pficm which is more metallic than many conventional metallic 
alloys. Thus varying x (and hence the electronic concentration) one can go from a good 
metallic state to an insulating state. The sodium tungsten bronze system (Na^WOa) has 
been well studied because of this composition driven metal-insulator transition seen in this 
systemic. This system has an additional feature that when W is substituted by Ta, the 
empty 5d orbitals of Ta compensates for the electrons being introduced by Na. This reduces 
the carrier density and also introduces a large disorder. One can thus see the effect of the 
disorder on the resistivity much like one does in a conventional alloy. On introduction of 
Tantalum in the tungsten site the chemical formula can be written as Na^Ta^/W^Oa. 

To summarize, the tungsten bronze system used by us has three distinct advantages : 
(1) the oxide is a simple oxide with a cubic structure, (2) the metallic state can be tuned 
to different electron densities by changing x and (3) the Ta substitution allows a way to 
introduce disorder into the system. The above three features thus gives us a good system 
in which we can make a comprehensive investigation. 
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III. EXPERIMENTAL DETAILS 



The samples studied were single crystals of composition x=0.9, 0.75, 0.63 and one sample 
with Ta substitution with composition x=0.6 and y=0.1 (x — y=0.5). The single crystals 
were grown by electrochemical growth from molten bath of sodium tungstate and tungstic 
oxide. Details of the sample preparation and characterization are given in the reference@0. 
We essentially followed the same technique. 

The effective carrier concentration n e j/ as a function of x or (x — y) is given by the 
relation n e ff = (x — y)/al, where a (the cubic cell dimension) follows the empirical relation 
ao=0.0820x+3.7845 A0 The values of n e ff calculated using the above relation and that 
obtained by Hall measurements^ are tabulated in table I. n e //(exp) is larger in comparision 
to n e jy(cal) by few tens of percent. In particular this discrepancy is larger for low n e yy. The 
reason for this discrepancy (while it is not very large) can be the assumption of the free 
electron formula for obtaining n e ff from the Hall coefficient. 

The resistivity of the samples were measured from 1.5K to 300K using four probe low 
frequency bridge technique^. The contacts to the samples were made by first depositing 
Au and leads were attached using silver paint or silver epoxy. Pre-cleaning of the surface 
before deposition of the Au is very crucial in this case. The dimension of the samples were 
less than 2mm x 2mm x 0.1mm. Because of the finite size of the contacts on the sample, 
the error in the measurement of the absolute resistivity is about 10%. The resistances of 
the samples measured varied between 1.0mfi-5.0fi. We could measure the resistance value 
with a precision of 5[iQ. 

The point contact spectroscopy was done at 4.2K in a dipper typeEl cryostat using Au 
as the tip to make the point contact. The tip was formed by electrochemical etching of Au 
wire using a solution of HC1, HNO3 and HF in the ratio 3:2:1. We used an ac bias of 5-10 
V for etching the Au wire. The I-V characteristics of the Au-sample point contact junction 
were taken by a ac modulation method. This directly gives the conductance (dl/dV) of the 
junction as a function of V. The bias used were upto 125mV. The junction was in contact 
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with the liquid Helium bath at 4.2K and hence there is no possibility of temperature drifts 
when the I-V characteristics were obtained. Typical junction resistances were adjusted to 
the order of 5-20 Ohms using a differential screw. The surface of the samples were carefully 
cleaned just before mounting the sample and were subjected to minimum exposure to the 
atmosphere to avoid contamination of the surface. 



IV. RESULTS 
A. Resistivity 

Figure 1 shows the resistivity of the samples as a function of temperature over the entire 
temperature range. The temperature dependence of the resistivity is like a metal and dp/dT 
> 0. It can be seen that the resistivity varies by a large amount by changing the value of 
x. In particular the Ta substitution can lead to a big jump in p. The value of the residual 
resistivity are given in table I. The residual resistivies of the samples agree well with the 
pervious reported valuesi'El. In figure 2 we show the variation of p(4.2K) as a function of 
x (for x > 0.5) obtained by other investigators^ on single crystals and data obtained from 
our single crystals. For samples with x<0.6 the residual resistivity generally varies within a 
factor of 2 depending on the extent of disorder. The exact value of p(4.2)K, however does 
not affect our experiement or the conclusions. 

The temperature dependence till 300K (i.e, Ap(300K)) for the low resistivity samples 
are nearly 15-25 pVt cm. While for the Ta substituted sample this can be as large as 88 fiQ 
cm. There is an interesting anomaly to the trend at x=0.75 so that the room temperature 
resistivity of this sample is actually lower than that of the x=0.9 sample. This arises because 
of partial ordering of Na atoms as has been seen previously^. We will see that this also 
affects Ap(300K). A detailed discussion on this issue is beyond the scope of the paper. 
Using a free electron model we estimate the mean free path at the two temperatures. The 
elastic mean free path Iq as well as the in-elastic mean free path /; n at 300K (giving rise 
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to Ap(300K)) are shown in the table I. The elastic mean free path varies strongly as x is 
varied and for the Ta-substituted sample l /a «2. For the most metallic sample (x=0.9), 
l /ao ~ 145. Thus in the range of the composition studied, the electrons are not on the 
verge of localisation. Most importantly we are in a position to assess the strength of inelastic 
scattering quantified by Z" 1 in the presence of varied l while keeping the system chemically 
unaltered. In comparision to variation of Iq with x the variation of U n is not that severe. 
Interestingly the Na ordering in x=0.75 leads to a drastic increase in l in . 



B. Point Contact Spectra 

In figure 3 as an example we show the second derivative d 2 I/dV 2 obtained from an 
experimental dl/dV-V curve for x=0.75. We also show the dl/dV curve for the sake of 
completeness. Here dl/dV decreases as — V — is increased and it is a sign of "metallic" 
point contact. The voltage axis directly gives the phonon energy in eV and as stated earlier 
d 2 I/dV 2 gives a measure of the phonon density of states. Essentially d 2 I/dV 2 arises from 
inelastic scattering of electrons by the phonons. For the other samples the curves were 
similar and we donot show them to avoid over crowding of data. As a comparison we show 
that data taken on Au (figure 3(a) inset) which acts as reference to show the differences in 
the PS of the two classes of metals. For conventional metals the maximum phonon energy is 
often < 20 meV. For oxides there is a large density of phonons at energies beyond 20 meV. 
It will be our attempt to see whether the phonons seen in the point contact spectroscopy can 
give rise to the necessary electron-phonon scattering so that the temperature dependence 
of the resistivity can be explained. This requires a qualitative analysis of the experimental 
PCS data to obtain the spectral function a 2 pc F(uo). It is important to note that the thermal 
smearing (^3keT=i meV at 4.2K) is much small compared to the phonon energies which 
are relevant for our work. As a result we will not apply any correction to our experimentally 
obtained spectra and treat them as essentially the zero temperature spectra. 
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V. ANALYSIS AND DISCUSSION 



A. Analysis of Point Contact Spectroscopy data 

The phonon spectrum can be obtained from the second derivative of the I-V character- 
istics of the point contacts from the relational 

y = -^n eff N(0MeV) (1) 

where T(eV) is the spectral function related to the strength of the inelastic interaction which 
gives rise to the features in d 2 I/dV 2 , fi e // is an effective volume of phonon generation within 
which the inelastic scattering takes place, iV(0) is the density of states at the fermilevel. 
When phonons cause the inelastic scattering F{eV) = a 2 c F{oj) which is related to the 
Eliashberg function a 2 F(u) by a form factori. In the measurement of this electron-phonon 
interaction by means of point contact, there is a considerable background which is due to 
the scattering due to non-equilibrium phonons. Taking this into account the final expression 
used for obtaining the a 2 c F(u) is given by 



const x ap C F{eV) + k 
V o 



dV 2 



pc 



de (2) 



In the case of a heterocontact the spectrum is given by the sum of phonon spectrum (PS) 
in both the electrodes^. We have used Au as one of the electrodes and the PS of Au has 
features between 10 and 17 meV (see inset of figure 3(a)). By subtracting the PS of Au 
one can easily identify the PS of the other electrode. The value of the constant (const) 
depends on the regime (ballistic or diffusive) we are working in and this is characterized 
by the Knudsen number K = lo/r, where r is the diameter of the contact region0. For 
hetero-junctions like ours where electrons in one electrode (like Au) have a very large mean 
free path, the Iq of the electrons in the electrode with higher p limits the Knudsen number 
K. In table II the relevant junction parameters are given. Thus the most metallic samples 
with K >1 (x=0.9 and 0.75) lie in the ballistic regime, the x=0.63 sample with K ~ 1 is in 
the transition regime and the Ta substituted sample with K <C 1 lies in the diffusive regime. 
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The difference between the two cases comes in the effective volume of the contact Fl e ff 
in which the inelastic scattering takes place and which contributes to the non-linearity in 
the I-V characteristics^. For x = 0.63 and x=0.Q;y=0.1 samples, the necessary corrections 
were applied to take into account that the point contacts were diffusive. In the case of the 
ballistic regime, the volume £l e ff — 8r 3 /3 and in the diffusive case fi e // = 7rr 2 / /4, where Iq 
is the elastic mean free path of the electron. All the relevant numbers of the point contact 
junctions are given in Table II. Eqn. ||] was used to obtain the a^ c F{uj). The constant k 
was determined using the fact that the observed spectrum coincides with the background 
function for eV > /iu; maa; Ei. uo max was chosen to be 125 meV. The above integral equation 
has been inverted to obtain the ai c F{uo) which will be the total of the contributions from 
both the electrodes. Since ai c F{u) is independent of the contact details, we have subtracted 
the PS of the Au from this data which contributes a small peak at u ~ lOmeV and 17meV 
as shown in the inset of figure 3(a). The resultant is the phonon spectra of the sample 
alone. In figure 4(a) we have plotted the phonon spectra obtained by inverting the d 2 I/dV 2 
data for x=0.9, 0.75, and 0.63. The curves have been shifted for the sake of clarity. It can 
be seen that there are distinct features at energies between 40-80 meV. We were able to 
resolve the spectrum into 6 major peaks by using peak analysis tools. For one spectrum 
the contributions are marked by dotted lines and peak positions by arrows. The frequencies 
where the peaks occur were essentially unchanged when we vary x, but there was a small 
but distinct increase in the spectral weights of the peaks. In figure 4(b) we have shown the 
spectral weights {a^ c F{uj)) of some of the peaks marked in figure 4(a) for x=0.9, 0.75 and 
0.63. The peaks 4 and 6 show a constant increase as a function of x but peaks 2 and 5 show 
a dip at x=0.75. This sudden decrease at £=0.75 leads to lower electron-phonon coupling 
constant and thus can be related to the low Ap value observed for this composition (see 
table I). 

Before making any mathematical analysis we would first like to establish that the tem- 
perature dependence of resistivity Ap(T) or the inelastic mean free path U n really arises from 
the phonons particularly from those in the spectra beyond 20meV. For most conventional 
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metals the strongest contribution of the phonon spectrum comes from phonons within u < 
20meV. The oxides thus differ from the conventional metals in the fact that a large density 
of vibrational states extend upto lOOmeV. The electron phonon coupling constant A is given 



From the experimental curve we obtain a 2 c F(u). We can connect a 2 F{oj) and a 2 c F(uj) by 
a factor G (a 2 c F(u) = Ga 2 F(u)). The factor G depends on the contact characteristics and 
it can be calculated knowing the mean free path l and the contact dimension For the 
ballastic regime the value is a constant G = 0.25. For the diffusive regime the value is given 
by the relation G = ((K-l) /2K) + (l/ AK 2 ) xln(l+2if). The values obtained for our samples 
are given in Table II. Using this corrected a 2 F(u), we can obtain the value of A. If the 
inelastic scattering of electrons arise from the phonons then l~ n (which is a measure of the 
scattering strength) should increase monotonically with A. In figure 5 we plot /^(SOOK) 
as a function of A. This clearly brings out that is indeed closely related to A. This 
confirms that the source of the inelastic scattering are the phonons in particular those with 
uo > 20meV which have a very large spectral weight in the range u ~ 40-80 meV. The value 
of A for weak coupling metals like Cu, Na etc. is ~ 0.15. For strong coupling Pb, A = 1.34. 
The oxide superconductors mentioned in section I have a A of about 1.0. Comparing the 
values of A obtained in our samples, we can say that the electron-phonon coupling strength 
is weak to intermediate (f» 0.2-0.4). In the next section we make a quantitative evaluation 
of p from the phonon data and compare it with the experiment. 



The temperature dependence of the resistivity Ap(T) can be calculated using the formula 
derived by Ziman, in terms of the a 2 r Y{u) (electron-phonon transport coupling function)! 



by! 




(3) 
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B. Analysis of resistivity using the PCS 
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where u p is the plasma frequency. The transport coupling function differs from the Eliash- 
berg function a 2 F{oj) only by a factor analogous to the relation between a 2 F{oj) and 
a 2 F{uj). Work by Allen and co-workers0 provides a basis for using a 2 F{uj) instead of 
a 2 r F{oj) in the above formula. We thus replace a 2 r F{uj) by a 2 F{uj) in eqn.|[ As mentioned 
in the previous section, a 2 F{ui) = (l/G)a 2 c F(u), hence using the value of G given in table 
II, we can evaluate the temperature dependence of the resistivity using eqn |j. The only 
unknown parameter in eqn |4| is the plasma frequency uj p which we use as a fit parameter. 

The resistivity obtained for the samples, were fitted using the above relation for the 
temperature dependent part. The fit are given in figure 6(a)-(d). The maximum fit error are 
given in table III. It is interesting to see that the data can be fitted with the experimental 
a 2 c F{uj) over the whole temperature range for all the samples. The maximum fit error 
< ±0.5%. The fit error was random and showed no systematic deviation either at low T or 
high T. The fit thus can be considered excellent because the only free parameter available 
is the plasma frequency us p which appears as a multiplicative factor only. The prefactor u p 
obtained in the above relation from the p data is listed in table III (tOp(exp)). The table also 
contains the values of uo p as observed experimentally by electron energy loss spectroscopy 
(LEELS) which we call w p (obs)B 

For all the samples the value of uj p obtained from the resistivity data u p (exp) (eqn |J) 
and that obtained from LEELS u;p(obs) agree to better than 20% -30%. We think that this 
is a very good agreement given the fact that they are obtained from two widely different 
techniques. There is, however, a systematic difference. co> p (exp) is always lower than o; p (obs). 
One reason for this could be that while u; p (exp) is obtained from the bulk measurements, 
o; p (obs) is obtained from the surface sensitive LEELS. In these materials the Na has a 
tendency to diffuse to the surface. This makes the surface somewhat Na rich compared to 
the bulk and thus the n e ff at the surface can be somewhat larger than that in the bulk. We 
feel this makes o> p (obs) systematically larger than u; p (exp). 

For the sake of completeness we would like to investigate whether instead of analysing 
p(T) using the complete phonon spectra we use a discrete frequency approach. In this case 
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we take a Bloch-Gruneisun type of electron-acoustic phonon interaction characterized by 
a Debye temperature and an electron-optic phonon contribution characterized by a 
temperature The resulting p is given by 

5 B °I T 5 

pi T) J + B^'CJVZT)- 1 (5) 

Results of such a fit is given in the table IV. A typical fit shown in figure 7 for x=0.9 sample. 
For comparison we have plotted the fit error obtained using eqns f| and |5] in figure 8. It is 
clear that eqn |5| can be used to fit the p vs T curve but the extent of agreement is poorer 
compared to that obtained from eqn |j. The value of 6d obtained is close to what one expects 
for these oxides. From the specific heat data@ on these materials the values of 9r> are ~ 300 
- 400K. The value of /cb#e ~ 55 - 80 meV. This corresponds very well to some of the peaks 
in the phonon spectrum (figure 4a). 

Though the values of 9d and 9e are reasonable, the resistivity can definitely be better 
explained by the complete phonon spectrum rather than assuming the above expression |5| 
to be valid. 

C. Effect of Ta substitution 

It can be seen from figure 1 that even a small (10%) substitution of W by Ta leads to a 
large change in p(T). Between the £=0.63 and x=0.6 y=0.1 samples, n e // differ by ~ 20%. 
But the change in p is by a factor of more than 6 and that in Ap(300K) by a factor of 
nearly 2. The empty Ta orbitals at random sites thus act as a strong source of disorder in 
the system leading to substantial electron scattering which brings down the elastic mean free 
path l to just twice the cubic unit cell dimension. The comparision of x=0.63 and x=0.6 
y=0.1 sample will thus be like comparing a "metal" and a "substitutional alloy". 

Our concern in this paper is the temperature dependent part Ap(300K) which changes 
by nearly a factor of 2 and the inelastic mean free path Z in (300K) which is also quite small 
(w 27 A) in the Ta substituted samples. We would like to see whether the changes to a 2 F[uj) 
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brought forward by Ta substitution can also explain the temperature dependent part. This 
is an important test of our work that the Ap(300K) can be quantitatively explained by the 
observed a 2 F(u>) even when there is large substitutional disorder and the elastic mean free 
path is very small. 

In figure 8 we show the a^ h F(uj) for the two samples. Ta has very similar mass as W. In 
the region u ~40 - 60 meV there is a large enhancement of a 2 F(uj) and the peaks have moved 
towards higher uo. The enhancement of the oc^ h F(u) in the region u ~ 40 - 60 meV makes 
the A go up by a factor of 1.5. This causes the U n to go down by a factor of 1.75. Thus we see 
that the resistivity enhancement Ap(T) can be explained by the enhancement of electron- 
phonon interaction in the Ta substituted sample. This can be put into a simple description 
that whenever a large change occurs in the phonon spectrum a comparable change is also 
seen in the inelastic mean free path. Our above extensive analysis quantitatively establishes 
the fact that the temperature dependence of p in these oxides can be explained by the 
electron-phonon interaction. These modes mostly occur in the range 40 - 100 meV. We next 
discuss the origin of these modes and establish that they are optic modes. 



D. The optic modes 

WO3 is a polar crystal which has a basis consisting of differently charged ions. When Na 
is added, the empty conduction band is filled but the polar basis lattice remains unchanged. 
NaaWOs can be considered as a polar metal, which we can define as polar substance (WO3) 
with a partially filled conduction band. The structure of the WO3 (Re03 structure), is a 
vacant structure, where the oxygen ions can oscillate with large amplitudes (lower energy). 
The frequency of this oscillations correspond to the optical modes of the lattice vibrations. In 
a polar crystal, the optical vibrations are accompanied by a displacement and a deformation 
potential leading to polarisation. The total polarisation leads to a displacement which causes 
a strong coupling between the electrons and the optical phononsS The optical mode in WO3 
corresponds to a temperature of 600K or 50 me vi. 

As Na is doped into the system, the Na 
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occupies the vacant space and the amplitude of the oscillation of the oxygen ions decreases. 
This increases the energy corresponding to the oscillations. This shifts the spectral weight 
of the different peaks towards higher values of uj (see figure 4a). The phonon frequencies of 
Na^WOs mostly arise from the W06 octahedra and thus are similar to WO3. The phonon 
dispersions in W0 3 has been calculated^ and we use this as an identification guide to our 
experimental spectra. In WO3 the highest calculated acoustic mode frequencies (at the zone 
boundary) is ~ 25 meV. The lowest peak in a 2 F(uj) occuring below 30 meV (i.e, peak #T) 
is likely to arise from the acoustic phonons. The lowest calculated optical mode frequency 
(at zone center) « 33 meV. Thus all the peaks of a 2 F(uj) for u > 30 meV are from optical 
phonons. At the zone center the calculated phonon frequencies are 33 meV and 77 meV for 
uj < 100 meV. At the zone boundary, the phonon frequencies are calculated to be around 
uj ~ 55, 75, 85, 115 meV where the exact position depends on the symmetry directions. 
The calculation of electron-phonon interaction parameters show that it is wave vector (q) 
dependent and also symmetry dependent^. The resulting a 2 F{uj) will have to be determined 
by an integration of all q using the phonon dispersion curve. Such a calculation is beyond 
the scope of this paper. We only note that in most symmetry directions the calculated 
electron-phonon interaction parameters reaches a high value for qir/a,Q > 0.5 - 0.75. As a 
result peaks in ct 2 F(uj) are expected to occur for values of phonon frequencies intermediate 
between those calculated for the zone center (q=0) and zone boundaries (q7r/ao=l). We find 
that the observed peaks in a 2 F(uj) indeed occur between these limits. Also the splitting into 
longitudinal and transverse optical modes leads to a maximum of 5 distinct frequencies in 
the range 30meV< uj <100meV. This also matches with our observation. This establishes 
clearly that the peaks in a 2 F(uj) indeed are related to the indentifiable optical phonon 
frequencies of the WO3 structure. In insulating WO3 the electron-optical phonon coupling 
is quite strong due to the polar nature of the material. However, on introduction of Na as x 
increases the material becomes conducting which leads to screening of the interaction which 
reduces the electron-optical phonon coupling. 

The Ta substitution in place of W leads to an increase in the a 2 F(uj) for uj lying in 
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the range 40-60 meV. The a 2 F(uj) remains unaltered in the other frequency range. This is 
rather interesting and allows us to draw some interesting conclusions. The phonon dispersion 
calculations show that the lowest optical phonon branches describe a displacement field in 
which the W(or Ta) atom moves in one direction and the oxygen atoms move in the other 
direction. This leads to a large dipole moment and large electron-lattice coupling. This will 
therefore depend strongly on the charge on the W-site and the screening associated with the 
atomic charges. When Ta is substituted for W, there is a hole on the W-site. This reduces 
significantly the screening and the coupling constant is strongly enhanced. This explains the 
increase in the value of the A in the case of Ta substituted sample (x—0.6 y—0.1) compared 
to the £=0.63 sample. 



VI. CONCLUSION 



In this paper we have brought out the role of optical phonon modes of simple perovskite 
oxides in the resistivity. The phonon spectra, better still the electron-phonon spectral fre- 
quencies were directly measured by point contact spectroscopy. We find that, there are 
distinct phonon modes at enegies between 40-80 meV in the phonon spectrum (a 2 F(u)) ob- 
tained using point contact spectroscopy. Using this phonon spectrum we are able to explain 
the temperature dependence of resistivity of the metallic sodium tungsten bronzes. We have 
obtained the values of A the electron-phonon coupling constant. Comparing this with the 
values obtained for other conventional metals and superconductors, we can conclude that 
the electron-phonon interaction in this system is in the weak coupling limit. It may not 
induce superconductivity but this has significant influence on the temperature dependence 
of resistivity. The paper is quantitative and clearly shows that the large temperature de- 
pendence of p seen in these class of oxides essentially arises from interaction with optical 
phonons. 
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FIGURE CAPTIONS 

Figure 1 Resistivity of the Na^TayWi^Os samples. 

Figure 2 p^K values obtained from our investigation and Ellerbeck et.al.0 
Figure 3 Point contact spectra. (a)d 2 I/dV 2 for x=0.75 obtained at T=4.2K. Inset shows 
the Phonon spectra of Au. (b) dl/dV obtained by modulation technique for the same sample 
(Ro « 50). 

Figure 4 (a) The point contact electron-phonon interaction function ai c F(uj) obtained 
using eqn |2|. The peaks have been marked following the multiple peaks analysis, (b) The 
spectral intensity of some of the peaks as a function of composition x. 
Figure 5 The inverse of inelastic mean free path at 300K as a function of the electron 
phonon coupling strength (A). 

Figure 6(a-d) The fit to the resistivity using eqn.[|for the samples. The only fit parameter 
is Up tabulated in table II. 

Figure 7 (a)Fit to the resistivity using eqn.[5] for £=0.9 sample, (b) Fit error (%) for the 
two fitting schemes ((a) eqn|| and (b)eqn.§ for x=0.9 sample. 

Figure 8 a 2 F(u) for x=0.63 and x=0.6;y=0.1 samples. The difference in the spectral 
intensities clearly indicate greater electron-phonon interaction strength (characterized by A) 
for the Ta substituted sample. 
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TABLE I 



Composition 


n eff (calc) 


n eff (exp)@ 


Po 


Ap(300K) 


lo 


U300K) 


X 


y 


/cm 3 


/cm 3 


pfl cm 


pfl cm 


A 


A 


0.9 




1.57xl0 22 


1.95xl0 22 


3.20 


25.3 


550 


70 


0.75 




1.32xl0 22 


1.60xl0 22 


6.95 


15.3 


288 


130 


0.63 




1.12xl0 22 


1.49 xlO 22 


47.01 


42.0 


45 


47 


0.6 


0.1 


8.91xl0 21 


1.18xl0 22 


330.10 


88.0 


7.5 


27 



TABLE II 



Composition 


Rt 




G* 




X 


y 


Q 






A 


0.9 




5.0 


5.0 


0.25 


110 


0.75 




3.0 


1.7 


0.25 


134 


0.63 




15.0 


0.5 


0.19 


90 


0.6 


0.1 


18 


0.09 


0.053 


83 



f : Junction resistance 
| : Knudsen number 

* : Correction evaluated using expressions from ref [6] 

b : Contact dimension calculated using r = {Apl/^Rf/ 2 for Ballastic regime 
and Wexler's interpolation formula for the diffusive regime, (p and 1 are of the 
sample, since the contact resistance is limited by the sample because of it higher 
resistivity) . 
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TABLE III 



Composition 


^p(exp) 


LU p (0hs)B 


Fit error 




eV 


eV 




x=0.90 


2.58 


3.14 


± 0.4% 


x=0.75 


2.83 


3.02 


± 0.4% 


£=0.63 


2.23 


2.97 


± 0.5% 


x=0.6 y=0.1 


2.13 


2.77 


± 0.3% 



TABLE IV 



Composition 


Po 


A 


B 


D 


e E 


Fit error 




fifl cm 


fifl cm 


fiQ cm 


K 


K 




£=0.90 


3.17 


74.55 


16.94 


356.4 


916.6 


± 4.0% 


x=0.75 


7.00 


42.37 


7.09 


334.6 


756.1 


± 2.5% 


x=0.63 


47.05 


67.21 


18.62 


252.4 


645.7 


± 2.5% 


a;=0.6 y=0.1 


330.14 


296.19 


33.94 


328.2 


855.2 


± 1.5% 
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TABLE CAPTIONS 

Table I The important parameters obtained from the transport measurements. 

Table II Parameters from the point contact junctions for all the samples used in calcu- 
lating the true a 2 F{uj). 

Table III Comparision of plasma frequencies obtained by our method and LEELS0 
along with the maximum fit error obtained using eqn.^ to fit the resistivity. 

Table IV Parameters obtained from the fit of resistivity to eqn || 
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